Thermodynamics of the planar Hubbard model 
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The thermodynamic properties: specific heat, entropy, spin susceptibility \s and charge susceptibility Xc 
are studied as a function of temperature and doping within the two-dimensional Hubbard model with various 
U/t = 4 — 12. Quantities are calculated using the finite-temperature Lanczos method with additional phase- 
averaging for a system of 4 x 4 sites. Results show that the entropy at low T reaches a maximum near half-filling 
at the electron density n ~ 1 ± 0.15 in the whole regime of studied U/t. The pseudogap in Xs{T) becomes 
clearly pronounced for U/t > 8 while \c shows a maximum close to half-filling. The relation of results to 
those within the t-J model and to experiments is discussed. 

PACS numbers: 71.27,+a, 75.20.-g, 74.72.-h 
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The Hubbard model is the simplest prototype Hamiltonian 
for correlated electrons. It has been and still remains the sub- 
ject of numerous theoretical investigations in connection with 
the metal-insulator transition the interplay between the 
magnetism and the itinerant character of electrons, and pos- 
sible superconductivity emerging solely from the electronic 
mechanism. A particular attention has been devoted to the 
two-dimensional model (2D) on a square lattice, expected to 
capture the physics of superconducting cuprates. A lot of ef- 
fort has been put into the numerical studies of the ground state 
properties, using various quantum Monte Carlo (QMC) meth- 
ods [|. 

On the other hand, there are rather few studies of the 2D 
Hubbard model at finite T > 0, in particular away but close 
to the half filling, i.e. at the electron densities n ~ 1. In the 
latter regime the minus-sign problem prevents the application 
of the QMC method at low T in large systems [0]. Gross fea- 
tures of the specific heat CV(T) have been obtained via the 
internal energy E(T) using the QMC [B]. Results reveal the 
evidence of at least two energy scales at large U/t 3> 1, the 
larger one representing the upper Hubbard band. The behavior 
at low T shows a marked difference between an insulator at 
half filling n = 1 with CV(T) oc T 2 , and an anomalous metal 
at finite hole doping rih — 1 — n > (or analogous electron 
doping). Within the metallic regime the QMC method was so 
far not able to reach temperatures below the exchange scale 
J ~ 4t 2 /U, which sets up an characteristic energy of spin dy- 
namics and is thus essential for establishing the low-T physics 
at low doping. The uniform spin susceptibility Xs(T) has also 
been calculated [Q| using QMC, with even larger restrictions 
(smaller systems) at finite doping n ^ 1, and by Dynamical 
Cluster Approximation |BJ|. On the other hand, low-T prop- 
erties of the Hubbard model with U 3> t are believed to map 
well on the properties of the t-J model which is projected on 
the basis space without doubly occupied sites. Several static 
and dynamic properties of the planar t-J model have been re- 
cently calculated and followed well into the regime T < J 
using the finite temperature Lanczos method (FTLM) [|| ^]. 
Two most relevant conclusions on the thermodynamic prop- 
erties of the 2D t-J model [K|] are: a) normal-state entropy 



density s(T < J) is maximum at the 'optimum' hole dop- 
ing rih ~ where n* h ~ 0.15 at J/t = 0.3, b) a pseudo- 
gap temperature T* (rih), experimentally (among alternatives) 
defined with the maximum in the uniform spin susceptibil- 
ity x(T) [0 0, shows up also in the t-J model where 
T*(rih) decreases with doping and vanishes at the 'optimum' 
one, c) even at quite low T <C J and in the 'underdoped' 
regime rih < n *h some thermodynamic properties are close to 
the behavior of a semiconductor-like nondegenerate fermion 
gas 0. 

Our aim is to obtain thermodynamic results within the pla- 
nar Hubbard model, which is numerically (for an exact diago- 
nalization approach) clearly more demanding relative to the t- 
J model. We list some relevant questions which we address in 
the following: a) are there any qualitative differences between 
the thermodynamic properties of the planar t-J model and the 
Hubbard model at large U /t, b) how does the entropy 'opti- 
mum' doping shift with decreasing U / 1, c) is there a pseudo- 
gap scale also at smaller U/t. 

We investigate the Hubbard model given by 



H 



-'£( 



C is C js 



He.) + U y]minii, 



(1) 



where c\ s (ci S ) and n is are creation (annihilation) and number 
operators for electrons, respectively, and the sum (ij) runs 
over pairs of nearest-neighbor sites. We limit our calcula- 
tions to U/t = 4, 8, 12, where values range from the modest 
U < W, smaller than the bandwidth W = 8t, to the strong 
correlation regime U > W. Note that the latter case cor- 
responds to the physics of cuprates where the spin exchange 
J ~ At 2 /U ~ 0.3 t. 

We study numerically the Hubbard model on a square lat- 
tice using the FTLM based on the Lanczos procedure 
of exact diagonalization and a random sampling over initial 
wavefunctions. The advantage in the case of thermodynamic 
quantities is that they can be expressed solely in terms of a 
grand-canonical average of conserved quantities (ks = 1), 
i.e., 

(/} = Tr/(/V e , S„ ) e -("-^)/T /Tte -(*- MJ v e )/T (2) 
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where N e , S z and /i refer to the number of electrons, the total 
spin and the chemical potential, respectively. In the case of 
quantities as in Eq.(Q), the FTLM does not require the stor- 
age of Lanczos eigenfunctions, but only of Lanczos eigenen- 
ergies e", where j = 0, • • • , M (M represents the number of 
Lanczos steps) while n = 1, • • • , R runs over random initial 
Lanczos wavefunctions. We refer for the details of the method 
to Refs.[[7| ^]. Using FTLM in the above way we are able to 
investigate the model on the lattice of JV = 4 x 4 = 16 sites 
with periodic boundary conditions. 

The main limitation to the validity of results comes from 
finite-size effects. The latter can be substantially reduced 
by employing the boundary condition (flux) averaging JT^]. 
In a system with periodic boundary conditions the latter is 
achieved by introducing the uniform vector potential 6 mod- 
ifying the hopping elements t — > tij — t exp(z# ■ fy). We 
use furtheron N t uniformly spaced phases 8 instead of a fixed 
9 = 0. In this way results are essentially improved at lower 
U < W. This is particularly evident for noninteracting elec- 
trons with U = 0, where results on small lattices otherwise re- 
veal pronounced finite-size effects. In this case, using Nt 3> 1 
most properties discussed here become exact even on a finite- 
size lattice. 

Still the main restriction in the thermodynamic validity of 
our results comes from finite-size effects which show up at 
T < Tf s where they start to dominate results [^J. In the par- 
ticular parameter space U/t — 4—12, the 'optimum' cases 
are at n ~ 1 ± n* h with n* h ~ 0.15 (coinciding with largest 
entropy s = s max ) where Tf s /t ~ 0.1 — 0.15. On the other 
hand, Tf s increases towards n = 1 and n — > 0, 2, respectively 
[R]. Since the properties of the Hubbard model ([[]) on a bipar- 
tite lattice are symmetric around half-filling we present results 
only for the hole-doped regime nh = 1 — n. > 0. 

Using Eq.(^) we directly evaluate within FTLM the elec- 
tron density n = (N e )/N, the entropy density s, expressed 
as 

s = ]xifl/N + ({H)-fi{N e ))/NT, (3) 

and the spin susceptibility \s = ({S z ) 2 ) /NT. Quantities cal- 
culated as functions of ^ and T can be consequently presented 
as well as in terms of n and T. Using above quantities we also 
evaluate the specific heat CV = T(ds/dT)^ and the charge 
susceptibility - electron compressibility Xc = (dn/d/j,)T- 

Let us first discuss FTLM results for an overall behavior 
of the specific heat Cy(T) (per unit cell), as shown in Fig. 1 
for U/t = — 12 in the whole relevant T regime. At high 
T > 0.5 t our FTLM results in general agree with those ob- 
tained previously with the QMC method [^J. The advantage 
of FTLM is that we can reach lower T ~ Tf s ~ 0.1 t, well 
below the exchange scale T <C J ~ At 2 /U. The main mes- 
sage of Fig. 1 is that Cy reveals the existence of (at least) 
two energy scales which are well separated for U ^> t, i.e. 
for U = 12 t. The upper maximum is related to excitations 
within the upper Hubbard band and is well pronounced near 
half-filling. For a larger doping, i.e. for n < 0.85, these exci- 
tations merge with the lower Hubbard band. At lower U — At, 



the upper maximum is only weakly present even at n = 1, 
and disappears at smallest available doping rih — 0.95. Note 
also that at U = 4 t, apart from n = 1, CV merges even 
quantitatively with the noninteracting result, U — (prop- 
erties at U = in Figs. 1-4 are calculated for an infinite 
lattice). When discussing the relation of presented results to 
those within the t-J model we point out that the upper scale 
(upper Hubbard band) is projected out in the latter so results 
for Cy are typically different for T > t. [§]. 
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Figure 1: Specific heat Cv (per unit cell) vs. T for various elec- 
tron densities n near half-filling and different U/t. U = result is 
calculated for an infinite lattice. 

In the following we focus on the lower energy scale which 
is essential for the understanding of quasiparticle and low-T 
properties. In Fig. 2 we show entropy density s as a function 
of electron density n for different U/t = — 12 as well as for 
few lowest T/t = 0.1 — 0.3. First observation is that U > 
leads to an increase of s, which is largest at an intermediate 
doping nh = nt ~ 0.15. As expected, results for U = 12 t 
are even quantitatively close to the ones within the t-J model 
[Q, with the corresponding J = 0.3 t where the maximum 
s has as well been observed at n* h ~ 0.15 and such a dop- 
ing has been identified as an 'optimum' one. We should note 
that such a characterization of 'optimality' does not seem to 
be in conflict with the usual one related to highest T c since 
experimentally in several cuprates the maximum in T c and in 
the entropy [|l3|] appear to be quite close in doping. Plausi- 
bly, n* h can be related to the most frustrated case where the 
kinetic energy of holes (preferring an ferromagnetic ordering) 
and the spin exchange (favoring antiferromagnetism) are com- 
peting and therefore one could expect n* h oc J / 1. Moreover, it 
is evident from Fig. 2 that the 'optimal' doping n* h ~ 0.15 is 



3 



quite insensitive to U in a broad range U/t = 4 — 12. 
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Figure 2: Entropy density s vs. electron density n for low T/t = 
0.1 - 0.3 and different U/t. 

In Fig. 3 we present results for the spin susceptibility Xa {T) 
for various dopings close to half-filling n = 0.8 — 1.0 and 
U/t = 0—12. We first note that here the phase averaging 
method brings substantial improvement. This is evident by 
comparing Fig. 3 with QMC results on the same 4x4 lattice 
[Q] obtained at a fixed phase 9 = 0. 

As expected, the onset of U > leads to an increase of 
Xs (T) at lower T < t. It is however more interesting to fol- 
low the development of pseudogap features with increasing 
U/t. One of experimental definitions of the (large) pseudo- 
gap temperature is related to the maximum Xs{T = T*) Q]. 
In fact, T* defined in this way matches well with other ex- 
perimentally established crossovers ^ |To[ ]. It has been found 
|0, [| that T* (rih) determined in this way within the t-J model 
matches well experiments. As foreseen from the mapping to 
the t-J model with J = 0.3 t, we find in Fig. 3 essentially the 
same behavior for the Hubbard model with U/t = 12. On the 
other hand, the pseudogap maximum becomes shallower for 
U/t = 8, although the location T*(n) does not seem to shift 
substantially. The pseudogap features disappear at U/t = 4. 

Let us finally comment on results for the charge suscepti- 
bility Xc — dn/dfi, as presented in Fig. 4. For noninteracting 
electrons at U = 0, Xc is essentially T-independent (except 
very close to n = 1, due to the van-Hove singularity) and 
is equal to the single - electron density of states at the Fermi 
energy Xc = Nf- Well away from half-filling, i.e. in the 
'overdoped' regime' n < 0.8, the effect of U > is only 
quantitative to reduce Xc- This can be attributed to an overall 
decrease of the effective density Af(e) due to the transfer of 
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Figure3: Spin susceptibility \at vs. Tforn = 0.8— 1.0 and different 

U/t. 



states into the upper Hubbard band. We also note in Fig. 4 
that at the same time U > leads to an even flatter variation 

of Xc(n). 

More indicative and challenging is the development within 
the 'underdoped' regime rih < with a pronounced T and 
doping dependence. Very close to half-filling n ~ 1, we 
are at T > Tt a dealing with chemical potential /j, within the 
charge (Mott-Hubbard) gap. A small density of charge carri- 
ers rih, <C 1 in this regime behaves as in a doped nondegen- 
erate semiconductor, (as established within the t-J model at 
low doping [[fl]]) where 



■(H - e„)/T 



(4) 



Consequently, we get Xc — n>h/T. Such a behavior is ev- 
ident in Fig. 4, quite universally for all U > and its va- 
lidity extends at lowest T up to nu ~ 0.1. A large increase 
in the maximum Xs at low T, being again rather insensitive 
to U/t, is a clear manifestation of strong correlations and of 
the increasing effective density of states Afp on approaching 
the metal-insulator transition. In fact, it has been claimed on 
the basis of the T = QMC results [|lj| that within the Hub- 
bard model the charge susceptibility diverges approaching the 
half-filling as Xc oc (1 — n)^ 1 / 2 . The latter is qualitatively 
consistent with the flattening of the chemical potential as a 
function of doping fi(n — ■> 1) in La2_ 2; Sr :I ;Cu04 observed via 
the ARPES measurements [|l4|]. Nevertheless, at given lowest 
T ~ Tf s ~ O.lt we cannot distinguish a scenario with an 
enhanced but finite Np at T = from a divergent behavior. 

Let us summarize some essential conclusions of the present 
study of thermodynamic properties of the planar Hubbard 
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Figure 4: Charge susceptibility \<J vs. n for T/t 
different U/t. 
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model: 

a) The FTLM seems to have advantages with respect to QMC 
and other numerical methods for the calculation of thermody- 
namic quantities away from half filling. The phase averaging 
method used in this study represents an essential improvement 
and to large extent reduces finite-size effects, in particular at 
moderate U < W. Using the FTLM and a phase averaging 
we reach in our study the low-T regime, i.e. T < J, where J 
is an effective scale where the spin exchange is fully active. 

b) At large U/t ~ 12 results for the thermodynamic quanti- 
ties match even quantitatively those of the corresponding t-J 
model (with J — At 2 /U) [0, [| in the low temperature T < t 
window. Excitations into the upper Hubbard band contribute 
significantly only at large T > t. On the other hand, for 
smaller U ~ 8 t both scales start to merge, and become in- 
separable for n 7^ 1. In this intermediate U regime there is 
still a qualitative but not a quantitative resemblance to low- 
T results within the t-J model. Finally, results for U ~ 4 t 
approach the behavior of noninteracting fermions. 

c) The effective exchange scale J seems to determine the 'op- 
timal' doping for the entropy maximum s(n = 1 — nf/) = max 
as well as the pseudogap scale T*{nh) in Xs(T*) = max. It is 
clear that only at large U we observe J ~ J ~ At 2 /U . On the 
other hand, for U < 12 t we see that the positions of extrema 
in s(n) and in Xs(T) are quite insensitive to U/t, indicating 
a rather constant J ~ 0.3 t as well as the 'optimum' doping 
nt ~ 0.15. The former fact can be understood in terms of 
less localized character of spin degrees, which leads to an ef- 
fective spin exchange interaction reduced relative to the large 
U expression J = At 2 /U . 



d) The pseudogap feature (maximum) in Xs (T) is well visible 
at U = 12 t, but remains only weakly pronounced at U/t = 8 
and finally vanishes for smaller U. This is consistent with the 
interpretation that the (large) pseudogap T* is related to an 
onset of short-range antiferromagnetic correlations, which are 
only weakly pronounced for U ~ 8 t away from half-filling 
n < 1. 

e) One expects also an analogous pseudogap in the specific- 
heat coefficient 7 (T) = C V (T)/T 0, where a depletion 
should appear at T < T* (n). It is evident that such an effect is 
present within the model, since near half-filling we see C v oc 
T 2 while at larger doping rih ~^ n* h we get C v oc T v with 
v < 1. 

f) We should note that the maximum in Xs (T) is not specific 
for the 2D Hubbard model, but seems to be generally present 
also in the ID model [|l5|]. Nevertheless, in a ID system there 
is no qualitative change in the character of low-energy (spin 
and charge) excitations on doping since the excitations have 
all the way a linear dispersion and consequently a nonvanish- 
ing 7(T — > 0), in contrast to a 2D system. 

g) Previous studies of thermodynamic quantities within the 
t-J model [Q] have shown that results (at J/t — 0.3) are 
even quantitatively in agreement with the experimental ones 
in hole-doped cuprates, in particular the doping dependence 
of the entropy s [jl3[ ], the spin susceptibility Xs l@ and chem- 
ical potential /i [[14|. Our results show essentially equivalence 
of the low-T behavior of the t-J model and Hubbard model 
with large U ^> t, hence the correspondence with experiments 
applies again. However, we have shown that many results do 
not change significantly in a broader range of U, i.e. there is 
even a quantitative similarity of s(t), nt, T*(nh) etc., so the 
agreement with experiments persists also in a broader range 
of U/t, 
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